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ABSTRACT 

We present an extension of the formalism recently proposed by Pepper & Gaudi to 
evaluate the yield of transit surveys in homogeneous stellar systems, incorporating 
the impact of correlated noise on transit time-scales on the detectability of transits, 
and simultaneously incorporating the magnitude limits imposed by the need for radial 
velocity follow-up of transit candidates. New expressions are derived for the different 
contributions to the noise budget on transit time-scales and the least-squares detection 
statistic for box-shaped transits, and their behaviour as a function of stellar mass is 
re-examined. Correlated noise that is constant with apparent stellar magnitude implies 
a steep decrease in detection probability at the high mass end which, when considered 
jointly with the radial velocity requirements, can severely limit the potential of oth- 
erwise promising surveys in star clusters. However, we find that small-aperture, wide 
field surveys may detect hot Ncptunes whose radial velocity signal can be measured 
with present-day instrumentation in very nearby (< 100 pc) clusters. 
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1 INTRODUCTION 

Open Clusters have long been used as laboratories to test our understanding of star formation and stellar evolution, as 
each contains a (sometimes large) sample of stars with relatively well-known and common properties (age, composition, 
environment) but spanning a wide range of masses. For the same reasons, since th e discovery of the first extra-solar planet 
around a Sun- like star in the field just over a decade ago (jMavor fc Queloall995l ). the possibility of discovery extra-solar 
planets in open clusters has been a tantalising goal. 

Open Clusters tend to be relatively distant, and their members relatively faint compared to the field stars usually targeted 
by radial velocity surveys. Most of the projects searching for extra-solar planets in Open Clusters therefore employ the transit 
technique, which is particularly well suited to dense stellar environments and has the additional advantage of providing 
a direct measurement of the planet to star radius r atio. Recent or ongoing transit surveys in Open Cluste rs include the 



USt APS (the University of St Andrews Planet Search. Istreet et aklliooi : ISramich et al.ll2005l: 



OC ( von Braun et allboosi ). PISCES (Planets in Stellar Cluste rs Extensive Searc h, iMocheiska et al 



Hood et al. 



200. 



20051). E XPLORE 

5l. bqoel) 



STEPSS 



(Surv ey for Transiting Extra-solar Planets in Stellar Systems, Burke et al. 20061 ) and the Monitor project ( Aigrain et al 
20071 ). 



Understanding the factors that affect the yield of such a survey is vital not only to maximise its detection rate, but also 
to enable the interpretation of the results of the survey, including in the case of non-detections , in term s of constraints on 
the incidence and parameter distributions of planetary companions. Recently, IPepper fc Gaudil l|2005a[ ) (hereafter PG05a) 
introduced an analytical formalism to estimate the rate of detection of exo-planets via the transit method in stellar systems. 
One particularly interesting result is the fact that the probability that transits of a given system in a given cluster are 
detectable, if they occur, is a very slowly varying function of stellar mass in the regime where the photometric performance is 
dominated by the source photon noise, but drops sharply with stellar mass in the background-dominated regime. This implies 
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that the number of detections to be expected from a given survey is roughly proportional to the number of stars with source 
photon counts above the sky background level in t hat survey. 

In a follow-up paper, IPepper fc Gaudil (|2005bl ) (hereafter PG05b) applied the aforementioned formalism to young open 
clusters, showing that transit surveys focusing on these systems have the potential to detect transiting Neptune and even 
Earth-sized planets, by making use of the fact that low-mass stars are relatively bright at early ages, and that their smaller 
radius gives rise to deeper transits for a given planet radius. This opens up the tantalising possibility of detecting transits 
of terrestrial planets from the ground, and what is more of doing so in youn g systems, where on e might obtain particularly 
interesting constraints on the formation and evolution of extra-solar planets (jAigrain et al.ll2007l ). 

The formalism of PG05a assumes that the photometric errors on each star and in each observation are independent of 
each other (i.e. that the noise is white). However, an o posteriori analysis of the detection thresho ld of the OGLE transit 
survey in the light of their RV follow-up observations of OGLE candidate transits (jPont et al.ll2005l ) demonstrated that the 
effective detection threshold is significantly higher than tha t expected for white noise only , suggesting that correlated noise on 
transit time-scales might be present in OGLE light curves. iPont. Zucker fc Queloj (12003) (hereafter PZQ06) since developed 
a set of methods for evaluating the amount of correlated noise on transit time-scales in the light curves of transit surveys, 
and applied them to the OGLE light curves to show the latter do indeed contain cor rela ted noise at the le vel o f a few mm ag. 
Similar analysis of light curves from other transit surveys (see e.g. Smith et al. 20061 and Irwin et al. 2007 . and Pont 2007 for 
an overview) has since shown that they are also affected by correlated noise at a similar level. Because correlated noise does 
not average out as more observations of a given transit event are obtained, as white noise does, it is generally correlated noise 
which dominates over white noise in determining the detectability of transits around all but the faintest stars in a given field 
survey. A number of effects can give rise to correlated noise, including seeing-dependent contamination of the flux measured for 
a given star by flux from neighbouring stars, pointing drifts combined with flat-flelding errors, and imperfect sky subtraction 
- some or all of which may be important in a given survey depending on the telescope/instrument combination used and the 
observing strategy. 

Photometry alone does not allow the mass of the companion causing the transits to be ascertained, and radial velocity 
measurements are thus generally needed to confirm the planetary nature of a transit event. As pointed out by PG05a, 
this effectively imposes an apparent magnitude limit for transit detections to be confirmable, as accurate radial velocity 
measurements of faint stars are extremely expensive in terms of large telescope time. 

As noise that is correlated on transit timescales reduces the detectability of transits around the brightest stars in a given 
survey, but the need to perform radial velocity follow-up implies that only around the brightest stars can transits be confirmed, 
both effects must be incorporated in the scaling laws used to estimate the number of detections expected from a given survey. 
The present paper attempts to do this by extending the formalism of PG05a to include correlated noise and by translating 
the magnitude limit imposed by radial velocity follow-up into a cluster-specific mass limit. 

Section [2] briefly sketches out the basics of the formalism of PG05a and describes how one or more additional noise terms 
representing correlated noise terms can be incorporated in it. The impact of these modifications on the noise budget on transit 
time-scales and on the transit detection probability as a function of mass are investigated in Section[3l Considerations external 
to the transit search itself, including radial velocity follow-up, are introduced in Section |4] Finally, the practical implications 
of the resulting formalism for Open Cluster transit searches are briefly explored in Section [5] 



2 INTRODUCING RED NOISE TERMS 
2.1 Overall formalism 

This formalism is described in detail in PG05a, and only its major characteristics are sketched out here, so as to allow the 
modifications implied by the presence of correlated noise to be made clear. 

PG05a compute the number of transiting planets with periods between P and P + dP and radii between r and r -\- dr 
that can be detected around stars with masses between M and M -f dM in a given stellar system as 



where A'dct is the number of detected transiting planets, A'^* is the total number of stars in the system, d^p/dr dP is the 
probability that a planet around a star in the system has a period between P and P + dP and a radius between r and r + dr, 
fp is the fraction of stars in the system with planets, Ptot(Af, P, r) is the probability that a planet of radius r and orbital 
period P will be detected around a star of mass Al, and dn/dM is the mass function of the stars in the system, normalised 
over the mass range corresp onding to A*. 

Following [Caudil (|200Gh . PG05a separate Ptot(Af, P, r) into three factors: 



Vtot{M, P, r)^Ptr{M, P)Vs/n(M, P, r)Vw(P) 



(2) 
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Vtr is the probability that a planet transits its parent star, Vs/n is the probability that, should a transit occur during a 
night of observing, it will yield a signal-to-noise ratio (S/N) that is higher than some threshold value, and Vw is the window 
function that describes the probability that more than one transit will occur during the observations. 
PGOSa's expression for the transit probability is used without modification 

P^^^^=(^^y\r'/^RP~y^ (3) 

where R is the star radius and a the orbital distance. 

The S/N of a set of transits is S/N = (Axtr)^''^, where Axtr is the difference in between a constant flux and a boxcar 
transit fit to the data. PZQ06 give the general expression: 

Ax?. = -2 = (4) 



EC., 

where d is the transit depth, is the uncertainty on the transit depth, n is the number of in-transit data points and C is 
the covariance matrix the in-transit fiux measurement^. If the noise is uncorrelated, the non-diagonal elements of C are zero, 
and ^ dj = cr| where at is the uncertainty on the i"' flux measurement. Additionally, if this uncertainty is constant, i.e. 
(Ji = (Jv, Equation Q further reduces to: 

Ax?. = r.( — (5) 

\0"w / 

which, for single transits, is equivalent to Equation (4) in PG05a. Note that the notation adopted here matches that of PZQ06, 
and thus differs that of PG05a,b. In particular, the symbols A'^tr and n, used in PG05a,b to represent the number of in-transit 
points and the number of transits respectivtely, are inverted here. We also use (Tq where PG05a,b used a, and d where they 
used 5. 



2.2 Modifying the detection statistic to account for red noise 

In an attempt to account for the saturation of the rms. noise level that is seen at the bright end of all transit surveys, PG05a 
introduced in their Section 4.3 thqe concept of a minimum observational error osya, which is added in quadrature to the error 
contribution (Jphot from the sky and source photon noise to give the error estimate aind for each data point: 

0"w = (o-phot + C^sys)'- (6) 

This expression for CTw is then simply inserted into Equation [5]) . 

However, detailed analysis of the light curves of various transit surveys fPZQOe. lPoii^boOTl l shows that they systematically 



contain noise that is correlated on transit timescales (2-3 h for a Hot Jupiter transit), i.e. the non-diagonal elements of the 
covariance matrix are non-zero. As a result, no longer decays as n~^^^ as expected for uncorrelated (white) noise. PZQ06 
propose a single parameter description of the covariance, assuming the noise can be separated into purely uncorrelated (white) 
and purely correlated (red) components, the former decaying as n^^^^ but the latter independent of n: 

2 

O-rf = h CTr (7) 

n 

where ctw and (Tr reresent the white and red noise components respectively. This single parameter description of the correlated 
noise assumes that the degree of correlation remains unchanged on all timescales up to the maximum transit duration. It is 
equivalent to approximating the covariance matrix with Cu = = + in the diagonal, dj = for two data points in 

the same transit, and dj = otherwise (see Section for the treatment of multiple trans its) . 

There is evidence that the correlation timescale in transit survey light curves is flnite ( Gould et al.ll20o3 ) . If this correlation 



timescale is shorter than the maxium transit duration, the above expression would underestimate the significance of long- 



duration transit events. This does not appear to be the case for light curves analysed by PZQ06 and iPond (|2007l ). where 
the noise remains correlated up to 3 h timescales. Nevertheless, it is interesting to investigate the impact of finite correlation 
timescales through a simple example. We consider a transit with a depth of 1% lasting 2h and observed with 15min time 
sampling, i.e. n = 8. In the white noise only case, if (Tw = 2mmag, ad = 0.71 mmag and Axtr ~ 200. If correlated noise 
is present, with (Tr = 1 mmag, the single parameter approximation gives ad ~ 1.22 mmag and Axtr = 67. If on the other 
hand the noise is correlated only over timescales up to 1 h or 4 data points, i.e. dj = for \i — j\ > 4, (Td = 1.10 mmag and 
Axtr ~ 82. In general, even if the characteristic correlation timescale of the noise is shorter than a transit duration, we expect 



^ One can show that the estimate of d which minimises the of the fit is the inverse-variance weighted average of the in-transit 
flux-measurements, and is thus the uncertainty on this average. 
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the single parameter correlated noise approximation adopted here to give an estimate of the transit significance that is much 
nearer to the true value than that obtained with the white noise approximation. 

We therefore adopt Equation (O for what follows. The white noise is assumed to be equal to the photon noise and 
modelled as of contributions from the source and the sky background: 



2 _ 2 ,2 _ 

^source ~r ^back 



(8) 



where and A'^b are the number of photons from the source and the sky detected in the photometric aperture. 

PZQ06 found that the distribution of the rms. of OGLE light curves over a typical transit time-scale of 2.5 h is consistent 
with a constant red noise level of (Jsys ~ 3 mmag, independent of apparent magnitude. Processing the light curves with a 
systematics removal algorithm such as Sys-Rem reduces Gsys to ~ 1 .5 mmag fo r the best objects. The work of the International 
Space Science Institute (ISSI) working group on transiting planets ( PontllioO? ) has shown that similar values are also typical of 
other surveys, with a correlated noise value of ~ 1.5 mmag for the best objects. We therefore adopt cTays ~ 1.5 mmag throughout 
the following calculations, which would correspond to very good ground-based photometry. While it is theoretically possible 
to reduce the level of correlated systematics further, this value is used because it is considered representative of the leading 
surveys currently in operation. 

In addition to this systematics term, a red noise component proportional to the white noise level (as a function of 
magnitude) is present in some surveys. This dominates over the systematics term in the domain where background photon 
noise dominates the white noise and is thus likely to be somehow associated with background subtraction. We therefore label 
it cTsuh- For the purposes of the present calculations, it is modelled as a term proportional to the background noise: 



0"sub 



fccTb 



(9) 



For the purposes of the present work we assume k = 0.2. This is the kind of values the ISSI team found for the correlated noise 



in the HAT and SuperWASP surveys. Most cluster surveys, su ch as the University of St Andrews Planet Search ( Street et al 



20031: 



Bramich et al 



( Mocheiska et al 



2005 



Hood et al. 



Hartman et al 



20051 : iMocheiska et al.ll2006l ) or Monitor (| Aigrain et all 12007 



2005h , EXPLQRE-OC Jvon Braun et allboosl). STEPSS iBurke et al.ll2006h. P ISCES 



Irwin et alj 



2003), have 



'better spatial sampling, and lower values of k might therefore be expected to apply, though preliminary analysis of test light 
curves indicates that k ~ 0.2 is also appropriate, if not an underestimate, for at least some of these surveys. In any case, this 
value is used here to illustrate the effects of noise of this type when it dominates the overall noise budget. 
The overall red noise budget is thus 

fc'iVb 



2_2,2 _2| 
(Jj. ^^sys ' ^sub ^sys ' 



(10) 



2.3 Multiple transits 

As correlated noise does not average out over transit time-scales, but does average out over repeated transit events, it is 
particularly important to consider the repeatability of transits in the detection process when one suspects correlated noise 
might dominate. In an appendix, PG05a derived an expression for Vs/n for multiple transits, which is based on the equation 

Atr(multiple transits) = A'tr Atr(single transits) = A^tr— (11) 

where A*'tr is the number of observed transits. This equation remains valid in the presence of correlated noise provided there is 
no correlation over long timescales (similar to the planet's orbital period). In PG05a, Pw has to be calculated separately for 
each value of the number of transits. This assumes that the number of data points in each observed transit is the same, and 
in practice one must therefore choose a minimum value for the number of data points in a partially observed transit above 
which that transit contributes to Vw, and below which it does not. 

PZQ06 provide a general formula which accounts for the number of data points in each observed transit: 

A?,(multiple transits) = n^t^^^Tv: — — - (12) 

where Wtot ~ '^f. ~ 1^'' is the total number of in-transit points and 

rif^y ni~ block 

is the noise integrated over the fc**^ observed transit. Vw then becomes a multi-dimensional quantity dependent on not only 
A'tr but for each A'tr, on the set of nt- As with the PG05a formalism, it must be evaluated numerically. 

In the present work, we make the assumption of homogeneous phase coverage, which allows us to ignore differences 
between Uk for the different transits, and enables us to (roughly) estimate the number of transits observed as a function of 
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period given the time sampling and survey duration. This can then be incorporated into Equation (jlip directly, therefore 
alleviating the need to compute Pw separately. One can approximate A^'tr as 

iV,. = = ^^^^ (14) 

where ttot is the total time spent on target, which is the product of the number of nights A^n and the average duration of 
a night tnight, and P is orbital period. Reality diverges strongly from the homoegeneous phase coverage assumption close to 
harmonics of the daily interruptions in the observations, but it follows the global 1/P trend (see Fig. 1 of PG05a). 



3 IMPACT ON THE NOISE BUDGET AND DETECTION STATISTIC 
3.1 Noise budget on transit time-scales 

Useful insights regarding the dominant noise sources, and how to mitigate those that have the largest impact on the tran- 
sit detection performance, can be gained by exploring the dependency of the the various noise components on the stellar 
parameters. We start from the following expressions, given by PG05a, for A's, A^b and n (recalling that n is called A^tr in 
PG05a): 
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) = v^l-&2 neq,o ( 
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(15) 
(16) 

(17) 



where M is the stellar mass; a is the index of the (power-law) mass-radius relation; fix is the index of the (power- 
law) mass-luminosity relation in the filter X under consideration; D is the telescope aperture; tcxp is the exposure time; 
d is the distance to the cluster; Ax is the extinction to the cluster. PG05a adopt the distance-dependent extinction law 
Ai = 0.5(d/kpc); Lx,g is the Sun's photon luminosity in the filter of interest, which we compute, following PG05a, as 

8n^cRlX~\AXx 

Lx,Q = , I TT=r\ 7i (18) 

exp (/ic/Ax.cfcTo) - 1 

where \x,c and AAx are the filter central wavelength and FWHM respectively; S'sky.x is the sky photon flux per unit solid 
angle; Q, is the effective area of the seeing disk, which we compute, following PG05a, as 
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'sec 



(19) 



where Ssee is the FWHM of the PSF; 6 is the impact parameter of the transit; 5t is the interval between consecutive measure- 
ments. In PG05a, 



(5t — tcxp -|- ircad 



(20) 



where tread is the readout time, which can be generalised to include any time spent off-target; A^s,o is the number of source 
photons in the aperture for a solar-mass star; Jieq,© is the number of points in an equatorial transit for a solar-mass star. 
Substituting for Ns and A'^b from Equations (|15|l and (|16|) into Equation (jSj gives 
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where we have introduced, following PG05a, 



C2 = 



47rd^S'sky,xfi 



JX,0 



which is the ratio of sky to source flux in the aperture for a solar mass star. Taking (Jr from Equation l|10p . a« 
Equation (|21|l and n from Equation (|17|) and substituting into Equation (O, 
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(21) 

(22) 
from 

(23) 



To simplify this expression we introduce two new constants 
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0.1 1.0 0.1 1.0 

M (Mq) M (Mq) 

Figure 1. Error budget on individual data points (left) and over a the duration of an equatorial transit (right) for the fixed and fiducial 
parameters of PG05a, assuming Csys = l.Smmag and k = 0.2. The black dashed, dot-dash, triple dot-dash and dotted lines show the 
source photon noise, background photon noise, background subtraction noise and systematics terms respectively, and the solid black line 
shows the total noise budget. The grey line on the right panel shows what the behaviour the total noise would have it all the components 
behaved as white noise. The grey vertical dotted lines show mark transitions between the different regimes, as defined in Equations I I27II 
to ll30t . 



Ci = A'a,0 rioq.Q 

which is the square of the background subtraction component for a Sun-like star, and 



s,© 



(24) 



(25) 



which is the total number of source photons collected during an equatorial transit for a Sun-like star. (Note that C3 is defined 
in PG05a but not used here.) Equation (|23p then becomes 

1/2 / , , ^ i-«-/3x 

(26) 
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I + C2 



M 



M_ 





The form of (Jcq, the depth uncertainty for equatorial transits {b = 0) and of the different terms that compose it is 
illustrated in Figure [T] (right panel). Also shown for comparison is the noise level per data point (left panel), or V(l). 

The relative importance of the red noise components is clearly enhanced over the transit time-scale. While the systematics 
term is the same for all stellar masses, the source photon noise is a steeply decreasing function of stellar mass, the background 
subtraction noise is even steeper, and the background photon noise is the steepest. There may thus be up to four noise 
regimes, starting with the systematics-limited regime at the highest masses, followed by the source noise-limited regime, the 
subtraction-limited regime, and finally the background noise-limited regime at the lowest masses. 

Equating each pair of components and solving for M yields the mass regimes in which each component dominates. This 
exercise was done by PG05a to obtain Afgky, the transition mass between the source and background noise- limited regimes, 
which for clarity we rename M^^h- 



Me. 







Given the two additional noise terms that have been introduced, the relevant transitions are now: 
2 



M,y,,, = id T^^^^ Mo 



(C4 Ci)^-'^«+--'^x Mq 
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(27) 

(28) 
(29) 

(30) 
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However, it is very easy for the source noise-limited regime to disappear altogether, because the source photon noise 
averages out over the duration of the transit whereas the systematics and background subtraction noise do not. Even if one 
ignores the background subtraction term, the source-limited regime disappears if M^^h > Msys,s- Given the set of fixed and 
fiducial parameters adopted by PG05a, the source limited regime exists only if CTsys < 0.5 mmag. Adopting a more realistic 
value of 1.5 mmag, there is a direct transition between the systematics- and subtraction-limited regime, which occurs at 

M,y,.,,b Mq. (31) 

\ fays / 

For the subtraction-dominated regime to exist requires k to be relatively large {k > 0.2). If this is not the case, there is 
a direct transition between the systematics and background limited regimes: 

M,y,,b=f^^j Me. (32) 



3.2 Detection probability ■ps/^ 



The detection probability is derived following the same method as PG05a, although we consider multiple, rather than single 
transits. 

A transit observed A'tr times is assumed to be detectable if it gives rise to a detection statistic Axtr ^ ^Xmin- K equatorial 
transits of a given system are detectable, one can derive a maximum impact parameter 6max up to which transits of such a 
system are also discoverable (PG05a). This arises because n — teqVl ^ b'^/St, where b is the impact parameter of the transit 
and teq the duration of an equatorial transit: 



/ 4P \^ 
\ttGm) 



(33) 



where R is the radius and Af the mass of the star. We have assumed that the planet radius r <^ R and ignored limb- darkening, 
which allows us to ignore grazing transits as both extremely rare and hard to detect, and to write 5 — (r/R)'^ where r is the 
planet radius. 

Therefore, the probability that transits of such a system are detectable, i.e. that Axtr ^ Axmim reduces to Vs/n = &max 
when integrated over b, assuming the impact parameters are uniformly distributed between and 1. However, the statement: 

Ax?. = AxcqVl-b^ (34) 

which is valid in PG05a, no longer holds here, because Axtr is no longer simply proportional to n. Instead, Equations (|4]) and 
Q imply: 

-1 / o \ -1 



iVtr 



(35) 



The expression for 6max is found by setting the left hand side of Equation (|35|l to Axmin ^-^^ solving for h. This yields: 



\ 
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A mm 



(36) 



Inserting the expressions for the different noise terms derived above gives: 



S/N 



-a-l3x 
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(37) 



This expression reduces, in the case of white noise only - i.e. when aaya and k vanish - to PG05a's Equation (15). A 



similar expression for Axcq, the detection statistic for equatorial transits, also ensues: 



AXcq = 



(«o) 



-/3x 



I + C5 



(a" ) 



+ fsys + C5 I -JT— 



-2/3 X 



(38) 



If Ax?q < Axmiii for a- particular star-planet system, transits of that systems are not detectable, whatever the inclination. 
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[t] 




Figure 2. Detection statistic Ax^q for an equatorial transit, for individual transits and assuming the systematics are white (as in PG05a, 
left) or for multiple transits and incorporating both red noise terms with Csys = 1.5mmag and k = 0.2 (this work, right). The line styles 
have the same meaning as in Figure [T] The grey horizontal dotted line marks the detection threshold Ax'^ ■ = 30 adopted by PG05a. 
The grey vertical dotted lines mark the lower and, if applicable, upper mass limits between which the transits are detectable. The grey 
vertical dashed linos mark the mass range where RV follow-up is feasible with FLAMES+UVES. 



The overall behaviour of Axoq as a function of mass is illustrated in Figure [2] (right panel) . Also shown for comparison 
is the single-transit A^eq obtained following PG05a, i.e. assuming the systematics are white (left panel). 

PG05a point out that, using a = 1 and jUi = 3.5, Axeq oc M^^^ and M^^^^ in the source and background noise-limited 
regimes respectively, which has the remarkable implication that the detectability of planetary transits is virtually independent 
of mass for all stars above sky, while it decreases rapidly for stars below sky. In white, source limited noise only, the number 
of detections from a given survey is thus roughly proportional to the number of unsaturated stars above sky. In the red noise- 
limited regimes, Axtr no longer depends on n, i.e. on h (provided the transit is not grazing, a given transit event contributes 
the same amount to the detectability no matter what the number of observations in that transit). Transits of a given system 
are thus detectable, whatever the inclination (i.e. the transit duration), if they are deep enough (5) and enough of them are 
observed (A'^tr). Using the same values of a and /?/ as in PG05a, A^oq oc and in the systematics and background- 

limited regimes respectively, with the remarkable implication that transits are detectable only around stars below a certain 
mass, determined by the systematics term. The reason is that the degree of correlation of the noise lowers the advantage of 
having longer transits and lower photon noise (a larger and brighter primary) compared to that of having deeper transits (a 
smaller primary). 

The combined effect of the different noise terms across the entire stellar mass range is to give rise to a peak in Axcq 
versus M, as illustrated by Figure (2] This immediately points to a potentially very simple way of evaluating wether a given 
type of planet is detectable at all in a given cluster with a given observational setup: find the 'peak mass', or stellar mass at 



which Axcq is maximised, by differentiating Equation (|38|) with respect to M and setting the derivative to zero 



C4 



M 

Mo 



-Px 



M 



+ C5(l + 2/3x) 



M 



(39) 



The solution of Equation (|39|l could then be plugged back into Equation (|38|l . to yield AXpcak- ^ AXpeak > ^Xmim detections 
are possible in the cluster under consideration. In that case, one can also compute the limits Miow and Mup of the stellar 
mass range over which detections are possible by setting the left hand side of Equation (|38|) to Axmin ^nd solving for M. In 
practice, both equations cannot be solved analytically in the general case. If one is interested in calculating precise values of 
■^Xpeakj Aflow and Mp, the simplest way is to to compute Ax?q for a range of M and find the quantities of interest numerically. 



However, as discussed in Section [3. II it is not uncommon for a single component to dominate over a significant portion of 
the mass regime. By considering dependence of Axcq on each of the components one at a time, one can obtain useful insights 
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into what limits the transit survey's performance, and what mass range will be accessible. Each of the source, background 
and background subtraction terms imply a minimum mass around which a given type of transit is detectable; 



Miow,b = {C[ C2) ^™r-T Mq 

where C( is a multiple-transit equivalent of PG05a's Ci 



n' — — 
- TF- - 



(40) 
(41) 
(42) 

(43) 



A^tr Nt^ Ns,Q ncq,Q y^e , 

and Miow,s and Miow,b are equivalent to PG05a's Mth.s and Mth,b- On the other hand, the systematics term implies a 
maximum mass: 



1 



(44) 



Note that PG05a's white systematic noise term also induces an upper mass limit detection, but it is typically larger than 
3 Mq. As long as A/up is above the peak of the mass function, correlated systematic noise will not significantly affect the total 
number of detections in a given survey. However, it does have the effect of preventing detections around the brightest stars, 
which are arguably the most interesting, because of their enhanced potential for follow-up. 

In general, more than one component contributes near the limits of the range of masses of stars around which transits of 
a given type of planet are detectable, and the expressions for Miow and Afup are rather complex. Again, an alternative is to 
compute each term in A^cq for an array of stellar masses and to find the values of M between which A^oq > Axmin- 



4 ADDITIONAL CONSIDERATIONS 
4.1 TurnofF mass 



Following PG05a, we compute the turnoff mass 



^ibol,0^ . 

where e is the efficiency of Hydrogen burning and (3 is the bolometric mass-luminosity index. 



(45) 



4.2 Saturation mass 

Also following PG05a's expression for the number of photons at the peak of the PSF of a given star, the saturation mass is 



l3x 



Lx,Q 



10 



OAAx 



(f) 



5sky,X(7pix 



cxp 



In 2 



9 

'sec 



(46) 



where A^pw is the full-well capacity of the detector and ^pix is the angular size of the pixels. 



4.3 Radial velocity follow up 

Radial velocity (RV) follow-up is necessary to confirm the planetary nature of any detected transits and to measure companion 
masses. In this section,we examine the range of stellar masses over which this is feassible for a planet of a given mass and 
period. 

PG05a used a fixed magnitude limit (V = 17 or V = 18) beyond which planets were considered undetectable by the 
radial velocity method. This is approximately suitable for planetary companions to Sun-like stars: it is extremely difficult to 
measure radial velocities with precisions of a few tens of m/s level beyond V ~ 18 even with the largest telescopes available at 
present ( Pont et al. 20051 ') . However, in cluster transit searches, many of the detections occur around lower-mass stars, where 



planetary companions may induce significantly larger radial velocity modulations, and a more detailed treatment is needed. 

For a star of a given magnitude, the minimum detectable radial velocity amplitude i^min is highly instrument depen- 
dent, and we examine two representative telescope / instrument combination: the UV- Visual Echelle Spectrograph (UVES) 
coupled to the Fibre Large Area Multi-Element Spectrograph (FLAMES) on the Very Large Telescope (VLT) - hereafter 
FLAMES+UVES - and the High Accuracy Radial velocity Planet Searcher on the 3.6 m telescope at La Silla - hereafter 
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HARPS. High precision measurements tend to be limited by instrument stability rather than by photon noise, in the sense 
that, if deemed interesting enough, a given (short-period) object can be observed as long as necessary, binning the phase-folded 
measurements to reduce the photon noise contribution. However, for each telescope / instrument combination there is also a 
magnitude limit Frv, beyond which the signal-to-noise ratio achievable in a single exposure drops below a critical level, and 
high precision measurements are no-longer feasible in reasonable exposure times. As the spectral region used typically covers 
the V and 7?-bands, Y should be either V or R, depending on which filter the object under consideration is brightest in. 
The radial velocity semi-amplitude induced by a given planet scales as 

K oimP-^^^M-'^'^ (47) 

where m is the planet mass and we have assumed that m <^ M and that the inclination of the system is edge-on. All planets 
giving rise to K > Jfmin are then assumed to be detectable around stars with apparent magnitude down to Frv, beyond 
which it is assumed that high precision radial velocity measurements are not feasible at all with the telescopes/instruments 
under consideration. For HARPS, we use F = V" and Yrv = 14, for FLAMES-fUVES we use Y = i? and Irv = 18. Both are 
relatively optimistic limits. This means that in each cluster, there is a lower mass limit 

MRV.min = 10"^-.O-^HV+51ogd-5+Ai,/2.5/3i- ^^g) 

where F is i? or V , below which no radial velocity measurements are feasible with a given instrument, and above which the 
minimum detectable planet mass is 

/ P y ( M V^' 

where mrd = MNoptuno for HARPS and Afjupitor for FLAMES -f-UVES. If considering a particular planet mass m across a 
range of stellar masses, one can derive a maximum stellar mass MRv,max around which such a planet produces a detectable 
RV signal by setting Afniin in Equation (|49p to m: 

/ „ N 3/2 / p \ -1/2 

MRV.max = Mo (50) 

V rrircf / \ 3 days J 

which is independent of the cluster properties and depends on the planet mass and period only. Planets with mass m and 
period P can thus be confirmed by radial velocity with present observational means only if they orbit stars with Mav.min < 
M < Mrv , max . 

Note that, for the sake of simplicity, we have ignored a number of important factors, including morphological differences 
in the spectra of stars of different types and the impact of rotation, which broadens the lines and degrades the radial velocity 
precision. 



5 APPLICATIONS 

One can roughly evaluate the mass range [Mmin; Mmax] within which planets of a given radius and period in a given cluster 
produce detectable transits and RV modulations: 

Mmin — max (A/low, AfRV,miii) (51) 
Afmax = min (Afup, Afto, AT^at, AfRV.max) (52) 



5.1 PGOSa's fiducial cluster 

Going back to the fiducial cluster of PG05a, under the relatively optimistic assumption that o-gya = 1.5 mmag, the detection of 
transits alone for a 1 Afjupiter planet in a 2.5-d orbit is possible around stars with masses between 0.28 and 1.49 Mq. However, 
using FLAMES-I-UVES on the VLT, for which we assume that A'min corresponds to a Jupiter mass planet in the same orbit 
around the same star and that -Rrv = 18, A/Rv,min = 1.13 Mq and AfRv.max = 1.22 Mq, so that the mass range where such a 
planet can be detect via transits and radial velocity is only 0.11 Mq. 

The combination of correlated systematics and follow-up requirements imposes very stringent limits on the potential of 
transit surveys in open clusters. In practice, it implies an even stronger dependence on cluster distance that illustrated in the 
bottom right panel of PG05a's Figure 8. 



5.2 Example galactic open clusters 

In a subsequent paper, PG05b applied the formalism of PG05a to a number of well-studied Galactic open clusters and, on this 
basis, made the prediction that close-in Neptune- or even Earth-sized planets should be detectable in some of these clusters 
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Name Distance Age Aperture tcxp Afio„ (Cause) Mup (Cause) AM 

(pc) (Myr) (m) (s) (Mq) (Mq) (Mq) 



Hyades 




46 


625 


1.8 


15 


0.15 


(lim) 


0.55 


(sat) 


0.47 


Praesepe 




175 


800 


1.8 


45 


0.15 


(lim) 


0.87 


(sat) 


0.79 


NGC 2682 


(M67) 


783 


4000 


3.6 


45 


0.15 


(lim) 


1.36 


(TO) 


1.28 


NGC 2168 


(M35) 


912 


180 


3.6 


45 


0.15 


(lim) 


1.49 


(sys) 


1.41 


NGC 2323 


(M50) 


1000 


130 


3.6 


45 


0.15 


(lim) 


1.49 


(sys) 


1.41 


NGC 2099 


(M37) 


1513 


580 


3.6 


45 


0.15 


(lim) 


1.49 


(sys) 


1.41 


NGC 6819 




2500 


2900 


6.5 


45 


0.15 


(lim) 


1.49 


(sys) 


1.41 


NGC 1245 




2850 


960 


6.5 


45 


0.15 


(lim) 


1.49 


(sys) 


1.39 


NGC 6791 




4800 


8000 


6.5 


45 


0.15 


(lim) 


1.08 


(TO) 


0.92 



Table 1. Masses ranges over which transits of Jupiter-sized planets in 2d orbits are detectable in selected Galactic open clusters, using 
the observational parameters of PG05b. Columns 7 and 9 give the primary cause of the upper and lower limits (sat: saturation; TO: 
turn-off; sys: systcmatics; lim: lower limit of mass range considered). 



Name Aperture texp Mi^^ (Cause) Mup (Cause) AM 

(m) (S) (A/q) (Mq) (Mq) 



Hyades 




0.1 


30 


0.15 


(lim) 


0.79 


(sat) 


0.70 


Praesepe 




0.1 


30 


0.18 


(RV) 


1.22 


(RV) 


1.04 


NGC 2682 


(M67) 


0.8 


15 


0.57 


(RV) 


1.22 


(RV) 


0.66 


NGC 2168 


(M35) 


1.3 


15 


0.64 


(RV) 


1.22 


(RV) 


0.59 


NGC 2323 


(M50) 


1.5 


15 


0.66 


(RV) 


1.22 


(RV) 


0.57 


NGC 2099 


(M37) 


2.7 


15 


0.75 


(RV) 


1.22 


(RV) 


0.48 


NGC 6819 




5.0 


15 


0.84 


(RV) 


1.22 


(RV) 


0.38 


NGC 1245 




6.5 


18 


0.93 


(RV) 


1.22 


(RV) 


0.30 


NGC 6791 




6.5 


65 


1.07 


(RV) 


1.08 


(TO) 


0.01 



Table 2. Masses ranges over which transits and radial velocity modulations of Jupiter-sized planets in 2 d orbits are detectable in selected 
Galactic open clusters, using SuperWASP for the Hyades and Praesepe and the observational parameters of PG05b for the other clusters, 
and using FLAMES+UVES for radial velocity follow-up. Columns 7 and 9 give the primary cause of the upper and lower limits (RV: 
radial velocity). 



via transit surveys from ground-based 2- to 6-m class telescopes. If so, transit surveys in open clusters might not only enable 
the detection of planets around well characterised stars of known age and metallicity, but may also lead to the first radius 
measurements for terrestrial planets. It is therefore interesting to investigate the detectability of Jupiter-sized and smaller 
planets in these clusters in the presence of red noise. 

We use the same test sample of 9 clusters (the Hyades, Praesepe, M67, M35, M50, M37, NGC 6819, NGC 1245 and 
NGC 6791) as PG05b, from which we take the cluster parameters (distance, age, extinction) and the observational parameters, 
which are similar to those of PG05a except that the night duration is tnight = 8h, the telescope apertures are 1.8 m (Pan- 
STARRS), 3.6m (CFHT) and 6.5m (MMT) depending on the cluster (as selected by PG05b), and the exposure time is 
icxp = 45 s for all clusters except the Hyades for which igxp = 15 s. 

Table [T] shows the range of stellar masses between which transits Jupiter-sized planets in 2 d orbits produce detectable 
transits. For all but the most distant cluster, transits are detectable right down to the minimum stellar mass considered 
(0.15Modot, well below the limit of 0.3 Mq adopted by PG5b), i.e. the addition of a correlated background subtraction term 
component does not affect the results. For all but the nearest clusters, the upper limit comes from the systematics term 
(A/up = 1.49 Msun, independent of the cluster and observational parameters). For the Hyades and Praesepe, the upper limit 
is saturation with the observational setup considered here, but this can be raised by using shorter exposure times and / or 
smaller telescopes, w hich have the added advantage of providing wider fields of view. For example, the SuperWASP project 
PoUacco et aD (|2006l ) uses multiple 11cm apertures and has 13.5 as pixels and an effective bandpass similar to R (Ssky ~). As 



it cycles between fields, S = 8min. With the standard exposure times of 30 s, we obtain Msat ~ 0.79 Mq for the Hyades and 
1.71 Mq for Praesepe. These represent a significant gain, and hereafter we adopt these observational parameters for these two 
clusters. Note that one could decrease the exposure time for the Hyades to increase Maat further, but this would confiict with 
the primary goal of SuperWASP, namely to search for transits around field stars. Overall, correlated noise does not strongly 
affect the detectability of transits of hot Jupiters in these clusters. 

We now incorporate the limits induced by radial velocity follow-up with FLAMES+UVES in the calculations. The results 
are shown in Table [2] Radial velocity or turnoff are now the limiting factors in almost all cases, and imply a stronger distance 
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Name Aperture foxp Mj*^^^ (Cause) Mj'^^ (Cause) A/up (Cause) AM 

(m) (s) (Mq) (Mq) (A/q) (Ms) 



Hyades 


0.1 


30 


0.08 


(lim) 


0.28 


(RV) 


0.50 


(sys) 


0.42 


Praesepe 


1.8 


15 


0.08 


(lim) 


0.63 


(RV) 


0.52 


(sys) 


0.44 


NGC 2682 (M67) 


1.8 


15 


0.18 


(back) 


1.06 


(RV) 


0.52 


(sys) 


0.34 


NGC 2168 (M35) 


3.6 


15 


0.16 


(back) 


1.19 


(RV) 


0.52 


(sys) 


0.36 


NGC 2323 (M50) 


3.6 


15 


0.18 


(back) 


1.25 


(RV) 


0.52 


(sys) 


0.34 


NGC 2099 (M37) 


3.6 


15 


0.32 


(back) 


1.49 


(RV) 


0.49 


(sys) 


0.17 


NGC 6819 


6.5 


15 


0.39 


(back) 


1.70 


(RV) 


0.46 


(sys) 


0.07 



Table 3. Masses ranges over which transits of Neptune-sized planets in 2d orbits are detectable (1) and confirmable with HARPS (2) 
in selected Galactic open clusters, using SuperWASP for the Hyades and the observational parameters of PG05b for the other clusters. 

dependence of the planet yield than transits. It is interesting to note, however, that confirmed detections of transiting hot 
Jupiters are possible down to very low stellar masses in the nearest clusters. 

We also investigate the detectability of Neptune-sized planets, using HARPS for radial velocity follow-up, still with a 
period of 2d. For such planets, the systematics term implies an upper mass limit of M^p = 0.53 independent of the 
observational set-up. We use the same observational setup as before except for Praesepe, where a greater photon-collecting 
capacity than SuperWASP's is needed to offset the smaller planet radius, so we revert to Pan-STARRS. The results are shown 
in Table [3] We have omitted NGC 1245 and NGC 6791 because transits of Neptune-sized planets are not detectable at all in 
these clusters. The systematics term severely limits the maximum stellar mass around which transits of hot Neptunes can 
be detected, while the need for radial velocity measurements limits the minimum mass around which they can be confirmed, 
and it is only in the Hyades that confirmed transiting Neptunes are expected to be detectable. We stress that these limits are 
relatively independent of theobservational setup. 

For hot Earths, the systematics term implies a very stringent upper limit of Mup = 0.13 Mq, and the formalism adopted 
here also precludes radial velocity confirmation around any stars in the clusters condidered (although it may be feasible to 
detect the radial velocity signal from a hot Earth around a bright star using HARPS by observing many repetitions of the 
orbit). 



6 CONCLUSIONS 

Simple modifications have been made to the formalism of PG05a to account for correlated noise and the need for RV follow-up. 
These should lead to more realistic estimates of the efficiency of Open Cluster transit surveys, while retaining the analytic 
nature of the original formalism, which affords useful insights into the behaviour of the detection probability as a function of 
mass. 

Two types of correlated noise were considered: systematics, which are constant with apparent stellar magnitude, and 
background subtraction noise, which scales with the background photon noise level. The latter behaves in a similar fashion to 
background photon noise itself, though its contribution to the total noise budget on transit time-scales has a slightly less steep 
dependence on the stellar mass, and therefore it does not significantly modify the yield of a survey unless extreme assumptions 
are adopted. However, the former implies a detection probability that steeply decreases with increasing mass and therefore 
curtails detections at the bright end. This effect is much stronger than the loss of sensitivity implied by a white minimum 
observational error of similar magnitude. 

In the course of evaluating the impact of correlated noise on the detectability of transits, we made a number of simplifying 
assumptions, and these should be borne in mind when comparing the predictions of the present formalism to the yield of 
real cluster transit surveys. First, we have assumed that the noise budget is the same for all stars of a given magnitude, 
and that every data point in a given light curve is affected by the same noise level. In fact, both white and correlated noise 
typically affect some objects and/or nights more than others, as they depend on factors which vary from object to object (e.g. 
crowding, position on the detector, colour) and time (e.g. weather, instrumental problems). Additionally, the way we compute 
the number of observed transits does not take into account the very strong features close to integer multiples of a day that 
are present in the window function of most ground-based surveys. As a result, while the scaling laws derived here apply for 
the majority of the objects in a given survey, the most significant detections in a real survey may well occur in special cases 
where the time sampling and the noise characteristics were particularly favourable. 

On the other hand, the radial velocity modulation induced by the companion in the primary, in order to measure the 
companion's mass, is only detectable given present day instrumentation over a certain stellar mass range which can be close 
to, if not above, the maximum mass implied by the systematics term for typical targets and observational set-ups. Thus, 
even though correlated systematics may not affect the yield of Open Cluster transit surveys significantly in terms of transit 
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detection alone (because transits usually remain detectable around stars close to the peak of the mass function) , it has a very 
serious impact on the yield in terms of transits whose planetary nature can be confirmed and where the companion mass 
can be measured. While the specific colour-magnitude relation followed by the members of a given cluster may enable one to 
exclude many of the astrophysical false positives which affect all transit surveys without actually detecting the radial velocity 
modulation of the primary, the scientific impact of any detection of a transiting planet will be significantly lowered if its mass 
cannot be measured. 

To illustrate a possible application of this modified formalism, it was applied to a selection of well-studied Galactic Open 
Clusters, which were used by PG05b to show that transits of Hot Neptunes, and even Hot Earths, should be detectable from 
the ground in nearby young Open Clusters. While correlated noise alone has little effect on the detectability of hot Jupiters 
in these clusters, we find that radial velocity follow-up severely limits the minimum mass around which their masses can 
be measured, which makes the confirmation of even Jupiter-mass planets in the more distant clusters difficult. Additionally, 
correlated systematics at the level of l.Smmag affecting all stars in a 20 night survey imply that transits of hot Neptunes 
are only detectable around stars with masses below 0.5 Mq. For such low stellar masses, the planetary radial velocity signal 
will only be measurable in very nearby clusters (< 100 pc) with present-day facilities. If hot Neptunes are abundant around 
M-stars, some could be detected by the combination of small aperture, wide field surveys such as SuperWASP and state of 
the art radial velocity instruments such as HARPS. 

The same level of systematics limits the detection of transits of Hot Earths to stars with masses below 0.13 Mq, irrespective 
of the cluster properties or observational setup. It is thus vital to achieve lower systematics (e.g. by going to space with CoRoT 
and Kepler) to detect transits of terrestrial planets, and particularly to detect them around stars bright enough that it may 
be possible to measure their radial velocity signal with future instrumentation. 

A general trend that emerges from this work is that the combination of correlated noise and RV follow-up requirements 
severely limits the choice of suitable target clusters, and effectively imposes a rather stringent distance limit. Additionally, we 
note that, for a given cluster, the optimal observational setup differs depending on the type of planet considered. 
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